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Abstract 

We study perturbation theory in certain quantum mechanics problems in which 
the perturbing potential diverges at some points, even though the energy eigenvalues 
are smooth functions of the coefficient of the potential. We discuss some of the 
unusual techniques which are required to obtain perturbative expansions of the 
energies in such cases. These include a point-splitting prescription for expansions 
around the Dirichlet (fermionic) limit of the <5-function potential, and performing a 
similarity transformation to a non-Hermitian potential in the Calogero-Sutherland 
model. As an application of the first technique, we study the ground state of the 
5- function Bose gas near the fermionic limit. 
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1 Introduction 



Perturbation theory is an old and often used method for studying problems in quantum 
mechanics which are close to some exactly solvable problems [|TJ. However, there are many 
models in which the traditional formulas of perturbation theory lead to divergent expres- 
sions for the energies even at fairly low orders, although the actual energy eigenvalues are 
known to be perfectly finite numbers. This situation typically (but not always) arises if 
the perturbing potential is singular, that is, if it blows up at certain points in the config- 
uration space. We are therefore forced to develop new kinds of perturbative techniques if 
we wish to expand such problems around the exactly solvable points. 

We will discuss two such techniques here. The first one is used in the presence of a 5- 
function potential in one dimension M , of which the ^-function Bose gas is a more general 
case ||. We will find that an expansion around the Neumann (bosonic) point, where 
the coefficient of the 5-function potential is zero, is rather straightforward. However, an 
expansion around the Dirichlet (fermionic) point, where the coefficient is infinity, diverges 
at second order in perturbation theory. We can obtain finite answers at that order by 
using a point-splitting prescription. After discussing all this in a one-particle model in 
Sees. 2 and 3, we generalise the discussion to the 5-function Bose gas in Sec. 4. 

The other technique discussed is used in the case where the perturbing potential has 
a scale-invariant inverse-square form; this leads to the wave functions vanishing as a 
power of the two-particle separation whenever a pair of particles approach each other. An 
important example is the Calogero-Sutherland model in one dimension j|, |5|, which has 
been studied extensively in recent years from many points of view (see || for a partial 
list of references). In Sec. 5, we show that this model can be studied perturbatively 
around the noninteracting bosonic limit by performing a similarity transformation of the 
Hamiltonian which takes into account the power-law form of the wave functions near 
coincident points. This leads to a perturbation which has a non-Hermitian form; one 
can then employ the usual Rayleigh-Schrodinger perturbation theory. As a byproduct, 
we prove that the energies get no contributions beyond the leading one or two orders in 
perturbation theory when the external (one-body) potential is either simple harmonic or 
zero; the wave functions however do get contributions to higher orders. 

Historically, the method of using a non-Hermitian perturbation was perhaps first used 
to study models of anyons in two dimensions J7| , in order to study the spectrum near the 
bosonic limit M. There is a vast literature on this subject, so we will not discuss it here. 
We should also briefly mention that there are other treatments of perturbation theory 
for singular potentials in two and higher dimensions, either in quantum mechanics || or 



in Chern-Simons field theories fL0 |. We will only consider one- dimensional models here 



which are technically simpler but have not been discussed as extensively. 

In this work, we will not be exhaustive in our list of examples of singular potentials. 
Further, certain parts of our discussion are not new, and all the models we will study are 
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exactly solvable. We should therefore explain the motivations behind our presentation. 
Our main aim is to introduce the reader to a class of problems in which the usual form 
of perturbation theory fails, and to indicate some ways of dealing with such problems. A 
more complete study of this subject would certainly be worthwhile but we will not attempt 
that here. Secondly, the models considered here would not remain exactly solvable if 
we introduced additional non-singular interactions or external potentials; however, such 
additional terms would not change the nature of the new perturbative techniques which 
are designed to address the singular parts of the interactions. It is therefore better to 
study the nature of these new techniques in exactly solvable models so that we can easily 
check that these methods do give the right answers. 

In our discussion, we will use the following expression from Rayleigh-Schrodinger per- 
turbation theory to quadratic order in a perturbation V. If H = Hq + V, and the exact 
normalised eigenstates ip n (0) and eigenvalues E n (0) of H are known, then the eigenstates 
of H to first order in V and the eigenvalues to second order in V are given by [|I[] 

E n = E n (0) + (^n(0)|^|^ n (0)> + 2^ ^77^ F77^ • I 1 ) 

l^n E n ([)) - El([)) 

If there is a degeneracy at energy E n (0), then we have to modify the procedure at first 
order in V by first diagonalising the matrix V\ n = (ipi(0)\V\ip n (0)) . However, we will not 
have to face this problem in any of the models discussed in this paper, because the matrix 
elements V\ n will always turn out to be zero for any two states for which E n (0) = Ei(0). 

In Sec. 5, we will also need the general form of perturbation theory at second and 
higher orders. The only information we will actually use is that, at all orders in V higher 
than the first, the change in the energy E n of the state ip n involves products of matrix 
elements like 

<VaO)|V|^(0)> (^(0)1^1^(0)) ...(^(0)1^(0)) (^(0)1^(0)) , (2) 
where k > 1, and all the intermediate states li, I2, h are different from n M. 



2 The (5-function potential in a box 

We begin with a particle which is restricted to move in one dimension in the range 
— 1 < x < 1. At the two end points |x| = 1, we demand that the wave functions should 
vanish. We impose a 5-function potential at the origin. The Hamiltonian is 

1 d 2 c 

H = " T- #2 + £ W > ( 3 ) 
2m ox z m 
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where c must be a real constant for H to be Hermitian. (We set Planck's constant h — 1). 
We will only consider the subspace of wave functions which are even functions of x, since 
the odd functions vanish at the origin and are therefore unaffected by the ^-function 
potential. The ^-function is equivalent to imposing the boundary condition 

^'(0+) = - V'(O-) = c0(O) , (4) 

where the primes denote derivatives with respect to x, and 0+ and 0— denote the limits 
x approaching the origin from the right and left respectively. (Note that c = gives us 
the Neumann boundary condition i/j'(0) = 0, while c = oo gives the Dirichlet boundary 
condition i/)(0) = 0). It is clear that this probem is exactly solvable; the (unnormalised) 
eigenstates and eigenvalues of H are given by 

ip n = sin [ a n (l - \x\) ] , 
a 2 

_ 

2m ' 

where ot n cot a n = — c . (5) 
The normalised eigenstates and energies for c = are given by 
ip n {x;0) = sin [ (n + — )7r(l - |x|) ] , 

E n (0) = 1 2m 2J , (6) 

where n = 0, 1, 2, .... If we expand the energies in (|J) around their values at c = 0, the 
leading order terms are 

E n (c) = -L [(n + i)V + 2c - / - +...]. (7) 
2m 2 (n + 2) 71 

For c = 00, the eigenstates and energies are 

ip n (x;oo) = sin [ (n + l)7r(l - \x\) ] , 
(n + 1) 2 tt 2 



EJoa) 



2m 



where n = 0, 1, 2, .... The energies in (||) can be expanded around their values at c = 00 
in powers of 1/c as 

(n + l) 2 n 2 2 3 

E n (c) = 1 ^ 1 -- + - + ... . (9) 

2m c c z 

We will now study how the expansions in (0) and (|9[) can be obtained by using per- 
turbation theory around c = and 00 respectively. Using the expressions in (|l|) with 
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V = (c/m)5(x), we recover the results in ([F[); at second order, we have to use the easily 
derivable identity 

4 1 

£ (n + \? - (1 + \) 2 = ~ (n + |) 2 ' (10) 

where the sum over / runs from to oo. 

We now have to discover how to perform perturbation theory around c = oo in powers 
of 1/c. We begin as follows. If ip n (x; c) and E n (c) are the exact results for the problem 
at any value of c, then the Feynman-Hellmann theorem says that with the perturbation 

V = (c/m)S(x), 

= I S ( x ) |^) = I |^(0; C )| 2 . (11) 

dc m m 

Using (|3J), this implies that 

</A '" " -|^(0+;c)| 2 . (12) 



d(l/c) m 

Therefore, the change in energy to first order in 1/c can be obtained by taking the per- 
turbation to be 

V = - 8"(x) . (13) 
2cm K 1 y ' 

Here, it is understood that when we compute (i[) n \V\i/) n ) , i.e., when calculating ip'(0;c), 
we must take the limit x — * either always from the right (0+) or always from the left 
(0— ). This is because ip{x] c) has a cusp at the origin for all c > 0; hence the two limiting 
values of ^'(0; c) are equal in magnitude but have opposite signs. 

Eq. ( |13D is one of the main results of this paper. It has not been discussed much 
in the literature before except for a brief mention in Ref. [Tl|] (there is a minor error 
in the relevant equation in that paper). It is straightforward to check that using this 
potential as the perturbation does give the correct first order change in energy if we use 
the c = oo results in (||) as the uperturbed wave functions. However, if use (|T3| ) in the 
second order expression for the energy in ([!]), we obtain a divergent sum. We therefore 
have to regularise the sum in some way. We do this as follows. From the first order result 
for ip n in ([!]), we find that for x > 0, 

ij)' n (x;c) = — (n + l)n cos [ (n + l)7r(l — x) } 

2(n + l)n + 1)2 cos [ (Z + 1)tt(1 - x) } 

(n + 1) 2 - (/ + 1) 2 ■ 1 ] 

It is clear that if we naively set x = in (13), the sum over I will diverge. 
We now use the result 

oo oo 

exp (ilx) = 2n 8 (x — 2np) , (15) 

l=— oo p=— oo 
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which follows from the completeness of the functions exp(ilx) in any interval in x of length 
27T. The right hand side of ( [T5D is zero in the neighbourhood of x = except exactly at 
the point x = 0. We therefore set the left hand side of the equation equal to zero in the 
limit x — > 0, although it diverges precisely at x = 0. This implies that 

3 

lim cos (/ + l)x = — - , (16) 

since I runs from to oo. We now subtract ( |16|) from the summation in (0), and then 
take the limit x — > 0+, to obtain a finite result 

V/(0+;c) = (-!)»(„+!).[! -? (§ + („ + !)* g__^_ 

= (-l)"(n +!)*[! - 1]. (17) 



Substituting this in (fT2[) , we obtain the correct second order result in 

Basically, we see that our regularisation of the energies at second order consists of 
using the first order wave functions not at the point x — 0, but rather at the point 
x = 0+ (or 0—). In the language of quantum field theory, this is called a point-splitting 
prescription. Note also that the second order change in the ground state energy (n = 0) 
in (H) is positive, although the formula in ([[]) would suggest that it should be negative 
(since E n (0) < Ei(0) for all I n); this peculiarity is a consequence of the regularisation. 

We will not pursue perturbation theory beyond second order near c = oo. It seems 
likely that a similar point-splitting prescription will also work at higher orders, but this 
requires a much more detailed analysis. 



3 Combination of (5-function and simple harmonic 
potentialss 

We now consider another exactly solvable example of the regularised perturbation theory 
described in the previous section. The motivation for studying a second example is to 
confirm that the point-splitting prescription is sound, and that it also works if the particle 
feels a potential in addition to the 5-function. We take the additional potential to be 
simple harmonic 0. Thus 

H = - J- ^ + - S(x) + \ mu 2 x 2 . (18) 
2m ox z m 2 

Once again, we restrict our attention to the subspace of even wave functions. The exact 
eigenstates and energies are given in Ref. ||. In particular, the energies are given by the 
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transcendental equation 



2^' 



s ^ 2 



From ([Tj]) , we can compute the first three terms in an expansion for E n around c = 
and c = oo. Near c = 0, 

£L 1 c (n — 1)!! c 2 /(n-l)!!\2 r n „ N , ,n 1 N1 ,„„ N 



2 y/nmuj nil 2^muj ^ n\\ ' 2 2 2 

where n = 0, 2, 4, we define 

p !! = p(p — 2)(p — 4). ..2 if p is even and >2, 
= p(p — 2)(p — 4)...l if p is odd and > 1 , 
= 1 if p = or - 1 , (21) 

and ^>(z) = d\ogT{z)/dz is the digamma function |12|] (not to be confused with the wave 
function ip n ). Near c = oo, 

where n = 1, 3, 5, .... The difference of digamma functions in fl20|) and (^) is given by 



V>(- + l)-^(- + -) = 2 log 2 if n = 



2 ' TV 2 2' 



k 



-1 

2 [ log 2 + 51 J ] if ^ is even and > 2 , 
k=i k 

n (_l)k 

2 [ - log 2 - ] if n is odd and > 1 . 

fc=i k 

(23) 



We will now derive the expressions ( p0|) and (|52 ) using perturbation theory. 
For c = 0, the exact eigenstates and energies are 

ij n = (mw) 1 / 4 H n (y) exp (-y 2 /2) , 

= ( n + - ) cu , (24) 

where n = 0,2,4, and H n is the Hermite polynomial of degree n in the dimensionless 
variable y = Xy/muj. To be explicit fT2| , 



/ n! y/2 ^ (-1)™ (2</) 2 ™ 

1 V 11 m=0 I 2 ~~ 



m)! (2m)! 



(25) 
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We now use the second-order formula in (|IJ) and discover that the results match fl2"0|) 
provided the following identities are true 



(J-l)!! 1 1 (n-1)!! r ,,n N , ,n 1 N , 

^nr^-^'*^ 11 -* 1 !^' 1 ' (26) 

where the sum over I runs over 0,2,4, .... For n = 0, this identity is given in Ref. JT3] , 
but for n > 2, we have not yet found it in the literature. However we have checked it 
numerically by computing the summation on the left hand side from I = to L = 2 x 10 7 



and setting the remainder (from I — L + 2 to I — oo) equal to y 2/ (ttL) by using Stirling's 
formula. We find that the two sides of (E^) match upto 1CT 9 for n = to 100. 
For c = oo, we have the exact results 

</v = (mw) 1 / 4 tf n (M) exp (-y 2 /2) , 

£ n = ( n + - ) u , (27) 
where n = 1, 3, 5, and H n is given by 

H n(y) = { t^tte ) E ^n^TW • (28) 



On using perturbation theory with (0), we recover the first order (1/c) result in (|22|) , 
but we again obtain a divergence at second order. We therefore use the point-splitting 
prescription and Eqs. ([!]) and (|12|) to obtain the second order energy 

uj v-+o+ c 2 ((^)!)% 3 / 4 2™ ^ 2'/2 {n -I) ' 1 J 

where we sum over I = 1,3,5,.... As before, we can extract a finite value for the limit 
x — > 0+ by using the completeness relation for Hermite polynomials 

E W #,.(!/) ^P [ - - ( y 2 + z 2 ) ] = 5 (y-z) . (30) 
p=0 / 

Setting z = 0, we get 



lim — — — ^-r. — 77- HJy) = . (31) 

2/^0+ „ (S ! 7T 1 /42P/2 PVy7 V ' 



p=0,2,4,... 



We use the differential relation H[(y) = VMHi^y) pf, and subtract (|TD from (g|). We 
then discover that (|29"D agrees with (^2|) provided again that the identities in (|26|) hold. 
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The results in this section indicate the general method of regularising second order 
perturbation theory near c = oo. Namely, we have to use the completeness relations which 
hold for the solutions at c = oo, as exemplified by the relations (p~5|) and (]30|) , in order to 
obtain a series whose sum is actually finite in the limit x — > 0+, although it diverges if we 
naively set x = 0. We then have to subtract this sum from the second order perturbation 
result in order to get a finite result in the limit x — > 0+. 



4 The (5-function Bose gas 

Having convinced ourselves that perturbation theory does work in the one-particle prob- 
lem with a 5-function potential, let us consider a many-body generalisation. We will study 
the ^-function Bose gas which has been studied extensively since it was first solved in Ref. 
using the Bethe ansatz. Although the exact results are available, it is instructive to 
see how they can be recovered by using perturbation theory near c = and oo. We will 
consider iV particles on a circle of circumference L. We will sometimes be interested in the 
thermodynamic limit N,L —>■ oo, keeping the density p = N/L fixed. The Hamiltonian is 

« = ~L E |? + ! E «<*.-*,)■ (32) 

Z '" L l<i<N UJj i l<i<j<N 

Since the particles are identical bosons, the wave functions must be completely symmetric. 
There are N\ possible orderings of the particle coordinates, given by < Xp x < xp 2 < 
■■■ < xp N < L, where (Pi,P 2 , ...,Pn) is some permutation of the numbers (1,2, ...,N). If 
the wave functions are known for any one ordering, say < x\ < X2 < ...xn < L, they 
are known for all other orderings by symmetry. 

Now it is clear that the model describes noninteracting bosons for c = 0. For c = oo, 
the wave functions vanish whenever any two particle coordinates coincide. We can then 
carry out the unitary transformation ipp —>■ (— l) p ipp, where ipp denotes the wave function 
for the ordering P, and (— l) p denotes the sign of the permutation P. Under this trans- 
formation, the wave function becomes completely antisymmetric, i.e., fermionic. Thus 
c = oo denotes a system of noninteracting fermions. Note that this unitary transforma- 
tion is only allowed if c = oo. At any other value of c, the symmetric wave functions 
do not vanish for Xi = Xj, and the transformation would produce antisymmetric wave 
functions which are discontinuous at those coincident points. 

For simplicity, we will limit our discussion only to the ground state wave function and 
the corresponding energy per particle. At c = 0, these are given by 

1 

= 



f = 0. (33) 
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We now consider the perturbative expansion around c = 0. At first order, the change in 
energy is found to be 

c(N-l) 



N L 

which is equal to cp in the thermodynamic limit. This agrees with the result in Ref. 
At second order, there are non-zero matrix elements between ipo and 



(34) 



■ipk-k = i 1 == ex P [ ik ( x i ~ x i) ] ' ( 35 ) 



(36) 



where k = 2-nn/L and n = 1,2,3,.... The second order formula in (]T|) then gives 

E<g ] _ c 2 (A-l) 
N 6 

This is indeed the correct answer (i.e., it agrees with the Bethe ansatz results) if we take 
the limit c — > holding A and L fixed. However, it is not the correct result if we hold 
c fixed and take the thermodynamic limit, as can be seen from the fact that the right 
hand side of ( |3~6"D diverges in that limit. The reason for this divergence is that the matrix 
element 

2cJN(N - 1) 

<^,-*IWo> = -^—j (37) 

can become comparable to or even larger than the energy denominator Eq — E^-k if 
we let L — > oo. Second order perturbation theory breaks down in that case. In the 
thermodynamic limit, a different perturbation theory exists due to Bogoliubov 0, and it 
gives the result 

%- = - ± P W . (38) 

A 671 

The presence of a fractional power of c clearly shows that this perturbation theory is 
rather different from the one we have been discussing so far. We emphasize that the 
divergence in our second order perturbation theory in the thermodynamic limit is not to 
be confused with the divergence which we have discussed in the Sees. 2 and 3; the two 
divergences are due to entirely different reasons. 

We will now consider perturbation theory around the free fermion point c = oo. Here 
we will again discover a divergence at second order perturbation theory which is similar 
to the ones discussed in the previous two sections, and the regularisation required to cure 
it is also very similar. However, unlike the expansion around c = 0, there is no additional 
divergence in the thermodynamic limit, and Eq/N does have an expansion in integer 
powers of the dimensionless parameter p/c. 

At c = oo, the ground state wave function ip in a particular ordering of the particle 



coordinates is given by l/vL N N\ times the determinant of the matrix 

Mij = exp [ ikiXj } , (39) 
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where the momenta fc, run over the N distinct values (vr/L) (1 — N, 3 — N, 5 — N, N — 1). 
(Note that for c = oo, the allowed values of the momenta are given by n2n/L if N is odd, 
and (n + \)2it /L if N is even ||). The ground state energy is given by 



E Q 1 ^ , o 7T 



2 



E = < " 1 ) ■ («) 



iV 2mN V 4 6mL 2 



Following the arguments leading upto (|T3"|), we can show that the perturbation around 
c = oo is described by the two-particle interaction 

V = - — Y 6"( x t - x, ) . (41) 

(The simplest way to derive this is to consider just two particles, and transform to the 
centre of mass and relative coordinates X = (xi + x<i)j1 and x = x\ — Xi respectively. 
The contact interaction then involves only the coordinate x which is a one-body problem 
as in the previous sections). Using this potential in first order perturbation theory yields 
the change in energy 



N 3mcL 3 



N ( N 2 - 1 ) . (42) 



We now proceed to second order perturbation theory. To simplify the calculation, we 
will consider the thermodynamic limit; thus we will replace J2ki by (L/2tt) J dk whenever 
convenient. It is useful to introduce the Fermi momentum 

k F = up , (43) 

so that all the momenta ki in the ground state lie within the range [—kp, kp}. We begin our 
analysis by observing that V has non-zero matrix elements between the ground state "0o 
and the excited states ipk',k'., where ipk'-,k'. differs from ipo in having exactly two momenta 
k[ and fc'- different from two momenta ki and kj] the other N — 2 momenta are identical 
in the two states. The translation invariance of the interaction ([4TD implies that we must 
have 

h + kj = K + k) , (44) 

where k\ and kj must lie outside the range [—kp, kp]. 
We then find that the second order energy is given by 

- 2 V ( k- - k ) 2 V ( k A ~ k 'i ) 2 (ak) 

N mc 2 NL 2 ,4 ^ ^ j k £ k , k 2 + k] - k? - kf ' ^ 45j 

which is again divergent since the sum over k[ or k'j runs over all values outside the range 
[—kp,kp]. Note that there are only three independent sums on the right hand side of 
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((45p since there is a constraint from momentum conservation ([44]); thus we can replace 
the sums by multiple integrals over ki, kj and kj of the form 

kp rkp roc 

dki / dkj / dk'j , (46) 

-kp Jki Jko 

where the lower limit of the k'- integral is given by 

k = kp if ki + kj < , 

= k F + ki + kj if ki + kj > . (47) 



We regularize the divergence in fl45j ) using the method described in Sec. 2. We will 
briefly indicate the procedure here without going through all the details. In the matrix 
element (i'o\V\4 , k' i ,k' j ) , let us consider one of the terms in V, say, 5"(x\ — x 2 ). After 
performing the integrals over £3, £4, xn, we are left with integrals over X = (xi + xz)/ 1 } 
and x = x\ — X2 oi the form J L dX/L times 

f 1 1 1 f 

/ ck(J(i) cos [-(k^—kAx] cos [-(&— A;,-) a;] = - / dx^i) [cos (k' j —k i )x + cos fe,-)a;] 
7 2 7 2'2j 7 J 

(48) 

appearing in the numerator of the last term on the right hand side of ([45]). If we naively 
perform the integral in (|5|) by setting x = 0, we obtain the constant 1; this leads to the 
divergence in (f|5[) when we integrate a constant over kj with the ranges shown in fl46|) 
and (0). To regularise that divergence, we first perform the integral over k'- for 1 ^ 0, 
and then take the limit x — > 0+. We use the result 

1 f°° , , 1 

- / dkj [ cos (kj — ki)x + cos (kj — kj)x } = — kp — - | ki + kj | (49) 

in the limit j^O. This follows from the fact that 

dk exp [ikx] =2 dk cos (kx) = (50) 
-00 Jo 



if x 7^ 0, and therefore 



/, x sin (i^x) 
dfc cos (kx) = i L 51 

K X 



which equals —K in the limit x — > 0+. 

After subtracting the integral in (^) from the right hand side of (f4"5|), we are left with 
the convergent expression 

fS 2 ^ 1 r k F r k F 



iV 47r 3 mc 2 p J-fc F Jfc, 



/ (ifci / c% ( fcj - kj ) I(ki,k 2 ) 

J —kp Jk; 



where I(ki,k 2 ) = 2fc^ + | ki + kj | — - (ki — kj) 2 f 

2 jfc 



fco J (Ay ki)(k'j fcj) 

(52) 
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On doing the integrals in fl5"2|) and using p3|) , we get 



p( 2 ) ^2 n 4 

- (53) 

Collecting all the results near c = oo, we find the first three terms in the energy per 
particle in the thermodynamic limit to be 

En tt 2 p 2 p p 2 , , 

-^ = -^1-2^ + 3^ + ... . 54 

N 6m c cr 

This agrees with the results in 0. 



5 The Calogero-Sutherland model 

We now turn to another technique for dealing with singular potentials. Consider the 
Calogero-Sutherland model in which particles move on a line (or circle), and interact 
pairwise through an inverse-square potential. This potential is so singular that the wave 
functions vanish whenever two particles i and j approach each other. This property can 
create difficulties for the usual kind of perturbation theory if we perturb around a model 
in which the inverse-square potential is absent and the wave functions do not vanish at 
coincident points. We will discuss two different versions of the inverse-square model, one 
with particles placed in a simple harmonic potential on a line, and the other with particles 
on a circle with no external potential. The modification of perturbation theory required 
in the two cases is similar, but it is convenient to discuss them separately. 



5.1 The Calogero model 

Consider the Hamiltonian for iV identical particles on a line Q 

H - _ J_ V— + A(A ~ 1} V - + — V x 2 (55) 

2m ^dx 2 + m (xi - x,) 2 + 2 ^ 

where A > 0. To make the problem completely well-defined, we have to specify that all 
wave functions vanish I whenever particles % and j approach each other. Since 

the singular two-particle interactions prevent particles from crossing each other, there is 
no way of comparing the phase of the wave function for one configuration with the phase 
for another configuration in which two particles have been exchanged. In other words, 
we have the freedom to work with any phase convention relating the different possible 
orderings of the particle coordinates. We will choose to work with completely symmetric 
wave functions. Note that the inverse-square interaction in (|55"D vanishes for both A = 
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and A = 1. However, the condition on the wave functions implies that A = corresponds 
to noninteracting bosons, while A = 1 corresponds to noninteracting fermions. The latter 
is true even though we are working with symmetric wave functions; the situation is very 
similar to the model in Sec. 4 with c = oo. 

It is particularly simple to solve the model at the two noninteracting values A = and 
1, and we can consider doing perturbation theory about either point. It turns out that 
perturbation theory around A = 1 is completely straightforward since all matrix elements 



of the perturbation (taken to be the inverse- square term in (55D) are finite. However, 
perturbation theory around A = is not so simple because even the first order result in 
(0) diverges. This can be seen from a simple counting of powers; if we define a relative 
coordinate x = x\ — X2, the wave functions for A = go to a non-zero constant as x — > 0. 
Then first order perturbation theory involves an integral like / dx/x 2 which diverges at 
x — 0. (The same power counting shows that perturbation theory around any value of 
A < 1/2 is divergent; the divergence is weakest (logarithmic) at A = 1/2 |T4|). 

We therefore proceed as follows for perturbation theory around A = 0. Let us denote 
the singular part of the wave functions at any value of A by 

A(A) = J] I Xi - Xj \ x . (56) 

i<j 

Instead of solving Hip = Eip, we will solve Hip = Eip, where ip> = A -1 -?/; does not vanish 
at coincident points, and 

H = A' 1 H A = H + V , 
1 ^_ d 2 muj 2 



where H ° = ~ 2^?^F + ^XX* 

V = - - £ (#-"#-)■ (57) 

m ^ Xi - xj OXi OXj 

We now have to derive the formulas for perturbation theory with the non-Hermitian 
interaction V given in (|57|). We can show that the required formulas are identical to 
those appearing in the usual Rayleigh-Schrodinger perturbation theory, such as Eqs. ([[]) 
and (0) 10. We only have to remember to use the expressions exactly as given in ([!]) 
and schematically as given in (0). For instance, we cannot replace the matrix element 
(ip n (0)\V\ipi(Q)) by the complex conjugate of (ipi(0)\V\ip n (0)) since V is not Hermitian. 
This will become clearer after the discussion below. 

The exact eigenstates and eigenvalues of the Hamiltonian Hq have the general form 
* n (xi ; A = 0) = P n (xi) G , 

G = ex P [ 7T 1^ ^i ] , 

z i 

4 0) = (n + j)u>, (58) 
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and P n is a symmetric polynomial in the X{ of degree n. In order to discuss perturbation 
theory, it is convenient to use the following result [TjJ. Any symmetric polynomial Q n (xi) 



of degree n can be written linear combination 

Q n ( Xi ) = J2 p i( x i) > ( 59 ) 
i 

where the symmetric polynomials Pi all have degree I < n and the wave functions PiG 
are eigenstates of H . 

Let us now consider first order perturbation theory. We note that 

V(P n G) = XN{N n - 1} P n G + Q t G, (60) 



where Qi is a polynomial of degree < n — 2. Using fl5"9|) , Q^G is seen to be a superposition 
of eigenstates of ifo with energy less than that of P n G. Therefore P n G is orthogonal to 
QiG. Hence the first order change in energy is just 

4- = • (61) 

One might worry about degeneracies at this stage. To see that this is not a problem, 
consider two polynomials P^ and P^ of the same degree n such that = P®G and 
?p( 2 ) = P^G are degenerate and orthogonal eigenstates of Ho. The orthogonality together 
with the form in fl60|) shows that both the matrix elements (ip^ \V\il)^) and (ip^ \V\ip^) 
are zero. It is therefore unnecessary to worry about degenerate perturbation theory at 
first order. We also have the general result that (ipi\V\tp n ) can be non-zero only if either 
ipi = ip n or if the degrees of the respective polynomials satisfy I < n — 2. This result will 
be used below. 

We now consider higher order perturbation theory. From the general form given in 
(0), it is easy to see that there are no contributions to the energies at any order higher 
than the first. This is because the individual matrix elements in (0) can only be non-zero 
if the degrees of the polynomials satisfy n < Ik — 2, Ik < h-i, h < h and li < n — 2. 
(These inequalities follow from the observation that the states ip^ may be identical to 
each other, but they all have to be different from ip n ). It is clear that all these inequalities 
cannot hold simultaneously; hence the product in (0) must necessarily be zero. 

Our perturbation theory has therefore yielded the well-known result that the energies 
only get a contribution at first order in A, and that this contribution is the same for all 
states as seen in (|61|). However the eigenstates can get contributions at any order in A. 

Before ending this subsection, we would like to emphasise that although the pertur- 
bation V is not Hermitian, its derivation from the Hamiltonian in ( |5^| ) guarantees that 
the perturbative changes in energies must be real. Note also that the eigenstates of H are 
orthogonal with the integration measure A 2 l\idxi, not with the measure Yiidxi. 
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5.2 The Sutherland model 



We now consider the inverse-square model on a circle of circumference L with periodic 
boundary conditions ||. The Hamiltonian for N particles is 

H= - — V — + V t . (62) 

2m i dx 2 m I? sin 2 [ j; { x i ~ x j ) ] 

We impose the same conditions on the wave functions at coincident points as in the 
previous subsection. To develop perturbation theory, we again perform a similarity trans- 
formation but now with a periodic function 

A(A) = I] I sin [ \ (xi - Xj ) } | A . (63) 

i<j 

This transforms the Hamiltonian into 

H = A' 1 H A = H + V , 

1 _ d 2 \ 2 ir 2 N(N 2 - 1) 

where H ° = ~ 2^ + — m 2 — ' 

and V = ■ - j £ cot [ j( Xl - x 3 ) ] ( A - A ) . (64) 
m L L oxi oxj 

Note that has a A dependent constant which will therefore appear in the energies of 
all the states. We now discuss perturbation theory with the interaction V. 

The (unnormalised) eigenstates and eigenvalues of the noninteracting Hamiltonian Hq 
are easy to derive. Consider a set of N integers {rii} = (m, n 2 , tin), each of which may 
be positive, negative or zero, and some of them may be equal to each other. Then we 
have 

ip{m} = exp [ i— {niXi + n 2 x 2 + ... + n N x^ ) \ + symmetnsation , 

Lj 

E in i} = U E n 2 . (65) 

The momentum of such a state is given by (2tt/L) J2i n i- 

In order to discuss perturbation theory, it is convenient to introduce the concept of 
squeezing ||. We say that a set of integers {n^} can be obtained from another set {rii} by 
squeezing if (a) two integers in the first set, say, n n! i and n'j, are different from two integers 
in the second set, say, n« and rij, but all the other N — 2 integers in the two sets are 
pairwise equal, and (b) v! i + n'j = rii + rij but — n'j\ < \rii — rij\. Thus the momentum 
of the states corresponding to the two sets are equal, but < £7{ ni } according to 
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We then find that the action of V on the states is given by 

2A 7T 2 

V ^ {ni } = ( — Y2 E I U i ~ n 3 I ) + E V>K> , (66) 

i<j {n^} 

where the summation over states runs over the various possible sets {?!■} which can be 
obtained from the set {rii} by squeezing; all these states have lower energy than ip{ ni } by 
the comment in the previous paragraph. We thus see that a matrix element (V'{n'-}l^lV'{r H }) 
can only be non-zero if {n'^} is either identical to {n^} or can be obtained from it by 
squeezing. Clearly, matrix elements between distinct but degenerate states are zero, so 
that we do not need to do any degenerate perturbation theory. 

The first order change in energy is therefore given by 

= £ I "< " | . (67) 

i<j 

Once again, there is no contribution to the energies from any higher orders in perturbation 
theory due to the product form in (0); at least one of the matrix elements in that product 
must be zero following an argument very similar to the one used in the previous section 
(with the concept of "degree of polynomial" replaced by "squeezing"). 

We have thus recovered the result that there are only two A dependent terms in the 
energy of any state. The term of order A is given in (|67| ) and it depends on the state, 
whereas the term of order A 2 given in (Bl) is the same for all states. 



6 Outlook 

There are several issues which could be addressed in the future. For the 5-function problem 
near c = oo in Sees. 2-4, we could study perturbation theory at cubic and higher orders 
to see if a similar point-splitting regularisation can be devised to obtain finite answers at 
all orders in 1/c. 

We can also study a generalisation of the 5-function problem to higher dimensions. 
Given a bounded region M with boundary dM, it is known |16| that the Hamiltonian 
H = —V 2 is self-adjoint (and therefore has real eigenvalues and orthonormal eigenvectors) 
if we impose the boundary condition on all wave functions, 

h-Vip = cip (68) 

where h is the unit vector pointing inward and normal to the boundary dM, and c can 
be any real number. Once again, we can study perturbation theory near c = (Neumann 
condition) or c = oo (Dirichlet condition). To do this, we first define a generalised 5- 
function and its derivatives which are non-zero only in an infinitesimal neighbourhood of 
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dM. We can then use expressions like (P or ( |T3"D where the perturbation V or V involves 
that ^-function or its second derivative. Once again, we obtain a divergence in second 
order perturbation theory near c = oo, and a regularisation is required to obtain finite 
answers. It would be interesting to study the precise form of the regularisation for an 
arbitrary region M. 

The method of non-Hermitian perturbation has already been used in two-dimensional 
models such as anyons |§ and systems with other types of two-body correlations []l7f . We 
can consider applications of this method to any model in which conventional perturbation 
theory diverges due to the singular nature of the wave functions in the vicinity of certain 
points in configuration space. We simply remove the singular part of the wave functions 
by a similarity transformation; we then obtain a new potential (generally non-Hermitian) 
whose matrix elements between non-singular wave functions need to be computed. The 
procedure will work if all such matrix elements can be shown to be finite. 
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